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The probability distribution in R, the proportion of neutral pions to all pions emitted in heavy 
ion collisions in certain kinematic regions, has been suggested as a signal of a disoriented chiral 
condensate (DxC). Here we note that the signature can be greatly enhanced by making suitable 
cuts in the data. In particular, we consider reducing the data set such that the k pions with lowest 
Pt are all neutral. We find that, given such cuts, (7?) can be substantially different from 1/3. For 
example, for a single DxC domain without contamination due to incoherently emitted pions, (R) is 
3/5 given the pion with lowest pr is neutral, and 5/7 given the two pions with lowest pr are both 
neutral, etc.. The effects of multi-domain DxC formation and noise due to incoherent pion emission 
^\ • can be systematically incorporated. Potential applications to experiments and their limitations are 

' briefly discussed. 
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I. INTRODUCTION 



Recently, there has been growing interest about the possibility of the formation of a disoriented chiral condensate 
(DxC) in heavy ion collisions Jl|-pT[. In an ultrarelativistic heavy ion collision, some region may thermalize at a 
^ ' temperature high enough so that chiral symmetry is restored in the region. As the system cools sufficiently rapidly 
back through the transition temperature, the chiral restored state is unstable as small fluctuations in any chiral 
direction (<j, 7?) will grow exponentially. This can create regions where the pion field has a macroscopic occupation 
number. It should be stressed that this scenario is not derivable directly from the underlying theory of QCD and 
contains a number of untested dynamical assumptions, principally that the cooling is rapid. While the failure of the 
system to form a DxC cannot be used to rule out that the system has reached the chiral restoration temperature (as 
the scenario described above is not derivable directly from QCD), observation of D%C formation would be a clear 
signal for chiral restoration at high temperature. 

There are clear signatures of D%C formation provided a single large domain is formed, containing a large number of 
pions. For example, one expects an excess in low px pion production as the characteristic momentum of a pion from 
a large region is small. Such a signal works even if multiple regions of D%C are formed provided each region is large, 
but it is not decisive since one could imagine some other collective low energy effects which produce low pt pions. 
On the other hand, since the pions formed in a DxC, being essentially classical, form a coherent state, this coherent 
state has some orientation in isospace, and all of the pions in the domain are essentially maximally aligned (given the 
constraints of quantum mechanics) and point in the same isospin direction. If there are a large number of pions in 
the domain, this implies a distinctive distribution of i?, the ratio of neutral to total pions in the domain |l|-^,|' 



/.<«) = ^ (i) 

In contrast, the distribution from uncorrelated emissions is narrowly peaked about 1/3 with a variance, (R 2 ) — (R) 2 = 
^ where J\f is the total number of pions, and approaches a delta function at R = 1/3 when AT — > oo. Since these two 
distributions are dramatically different, this provides a clear signature for single domain DxC formation, provided 
one can kinematically separate the pions from the DxC from other pions in the system. 

Unfortunately, this signature depends critically on the assumption that a single large domain of DxC is formed, 
which is a priori rather unlikely. One expects the domains are of characteristic size r, the exponential growth time of 
the pion fluctuation in the unstable chiral restored state |l|,[l4|] . Since the size of the fireball is much larger than typical 
QCD length scales, it seems unlikely that r would be of the size of the fireball, and thus formation of a single DxC 
domain is improbable. Formation of multiple domains of DxC, however, tends to wash out the R distribution described 
by Eq. (0). As the pions emerging from different domains cannot be distinguished kinematically, by the central limit 
theorem the R distribution will approach a normal distribution peaked at R = 1/3. This normal distribution may be 
distinguished from the normal distribution arising from uncorrelated emission; the case of multiple domains of DxC 
will have a substantially larger variance. 
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Unfortunately, there is an important practical limitation which makes it difficult to exploit the R distribution as a 
signature. Even under the most optimistic of scenarios, the total number of pions coming from D^C's will be a small 
fraction of the total number of pions. If one includes all pions produced in the reaction, the signal due to the pions 
from the D%C will presumably be overwhelmed. One can apply a low px cut to suppress the noise due to incoherently 
emitted pions. However, even with the low pt cut, the noise may still be severe, as both the signal and the noise 
peak at R = 1/3. In this paper, we suggest cuts which may dramatically enhance the signal-to- noise ratio. We study 
the conditional probability distribution of R given only for the events in which the k pions with the lowest pt are all 
neutral, and we will show that the expectation value of R is shifted away from 1/3. Since incoherent emission will 
result in a very narrow peak around R = 1/3, any such shifts should be easily observable. Moreover, one can make 
successive cuts by increasing the value of k, and enhance the signal in each successive step. 

This paper is organized as follows. We will start with the simple scenario and study single domain D^C formation 
in Sec. II. In Sec. Ill the effects of multi-domain formation and the noises due to incoherently emitted pions will be 
studied, while potential experimental application and limitations will be discussed in Sec. IV. 

II. DOMAIN WITH NON-ZERO ISOSPIN 

We are going to start with an unrealistic simple scenario and make it more realistic later on in our discussion. We 
will consider a single domain which is described by an isosinglet density matrix. (In general, a D%C is not a pure 
state; the wavefunction of the pion coherent state at the core of the "fireball" of the heavy ion collision is entangled 
with the energetic emission at the edge of the "fireball".) Such a state can be written as 

P= E c„/-^-|n,/,/ z )(n,/,/ z |, (2) 

n,I,I z ^ * 

with Y2 n i CnI = 1- Note that the coefficients c n i are real, positive and do not depend on I z (as the full state is 
assumed to be isoscalar) . The probability distribution of such a mixed state in the isospace is 

d 2 P{6,<f)) = ^c„//(2/+ 1) \ Y nA d ^)\ 2 sin6d9d4> = smdddd(f>/4n, dP{9) = sin 9d9/2, (3) 

n,I I? 

where the angles (9, 0) are defined such that a unit vector in isospace is (r x ,r y ,r ) — (sin 9 cos (p, sin 9 sin <p, cos 9) . The 
probability distribution is uniform, as the condensate is equally likely to point at any direction on the two-dimensional 
sphere, as demanded by isospin symmetry. 

The number operator of neutral pions in the condensate is given by 

n o = a l a o = ' a cos 2 9 = (n) cos 2 9, (4) 

with a — (a x ,a y ,a ) is a vector of hermitian annihilation operators which annihilate n x — (n + + 7T_)/v2, % = 
(tt + — 7r_)/\/2i, and t:q, respectively, and (n) — ^ n jC n in is the expected number of pions. Since R = no/{n) [], one 
can easily calculate the probability distribution of R. 

dP = \ sm0d6= id cos (9 = dVR. = \R~ l / 2 dR. (5) 

By defining dP = fo(R) dR (the subscript "0" stands for 1 = 0), one has, in the limit that (n) is large, 

MR) = ii?" 1 / 2 , (6) 

recovering Eq. (Q). The distribution is plotted in Fig. la. It is obvious that the shape is qualitatively different from 
the Poisson-Gaussian distribution due to incoherent emissions. The expectation value of R is 1/3, 



/ Rf (R)dR = l/3 
Jo 



(R) = / Rf (R)dR = 1/3, (7) 



1 Here we are assuming that the typical I is much less than (n) 
2 We are assuming that (n) 3> 1. 
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which has the simple interpretation that it is equally likely for the pion to be a ttq, 7r + or 7r_, and hence on average 
a third of the pions are neutral. 

This distribution is a consequence of the fact that we have assumed the D%C to be an isosinglet, a reasonable 
assumption on physical grounds. However, let's consider the distribution of R after the D%C emits a single ttq. The 
density matrix after the emission, which can be written as A aopaj, where A is a normalization constant, p is the 
density matrix defined in Eq. (|^) and ao annihilates a neutral pion. This new density matrix is not an isosinglet. It 
is straightforward to show that the probability distribution for this state is 



dP = 



\ sin 9 cos 2 f 



\R-R- l l 2 dR 



J \ sin 9 cos 2 6d6 fiiR. R-V*dR 



or equivalently, 



h(R) = f(R\ 



1st pion is neutral) = Rfo{R) j J Rf (R)dR = fi? 1 / 2 . 



(8) 



(9) 



The distribution ,fi(R) is plotted in Fig. lb, which is drastically different from fo(R)- The distribution is skewed 
towards the high end, while fo(R) is skewed towards the low end. Moreover, the expectation value of R is clearly 
pushed up: 



(R) 



i = [ Rfi(R)dR = 3/5. 
Jo 



(10) 



So we have arrived at the intriguing conclusion that, if the "first pion" emitted from a isosinglet D^C is neutral, 60% 
of the pions subsequently emitted from the D%C are neutral, a huge enhancement from the original expectation of 
33%. 




0.4 0.6 

Fig. Id 

FIG. 1. The probability distribution functions fk{R) for different values of k. Fig. la, b, c and d are fo(R), fi(R), fz{R) 
and fa{R), respectively. 
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This extraordinary statement certainly deserves more discussion. First, what is our criterion to decide which is the 
"first pion"? The answer is simple: it can be any criterion. It does not matter as long as it is a priori equally likely 
to be a 7To, a ir + or a 7r_. The derivation just depends on our removing a neutral pion from the isosinglet D^C. It 
can be the first pion emitted in time, or the last one emitted in time, or even the 17th emitted in time. The criterion 
can also be unrelated to the order of emission. For example, we can choose the "first pion" to be the one with the 
smallest polar angle. One can use any of these criteria to identify the "first pion" , and if it turns out to be neutral, 
then the R distribution of subsequent emissions is always given by f\(R), provided we are in the large number limit. 
However, this is only true in this idealized scenario, when all the pions are coming from a simple D%C domain. In 
reality, some of the pions come from incoherent emission, and if the "first pion" turns out to be incoherently emitted, 
the expectation value of R of the remaining pions is still going to be 1/3, not 3/5. As a result, we want to choose our 
criterion in such a way that the "first pion" is likely to originate from the D%C and not from incoherent emissions. 
Since D%C pions by hypothesis have low pr ~ l/£> where L is the size of the domain, a natural choice is to use the 
pion with the lowest pt as our "first pion" . 

After clarifying the meaning of the term "first pion" , we move on to discuss the physical origin of the modification 
of the probability distribution of R. In a nutshell, we are seeing the physics of (iso)spin alignment due to Bose con- 
densation. To illustrate the point, let us first consider the following apparently unrelated Stern-Gerlach experiment. 
Consider a large number of massive spin-1 particles, which for concreteness will be called deuterons. Initially they 
are all polarized along a randomly chosen direction n, which is a priori equally likely to be any direction in three 
dimensional space. In other words, S ■ n = for all the deuterons. Now let us pick one of these deuterons and pass it 
through a Stern-Gerlach spectrometer which measures S z , the spin along the z-axis. What is the probability that the 
measurement gives S z — 0? The answer is clearly 1/3, as the cases for S z = +1, and —1 are equally likely. On the 
other hand, if the measurement on the first deuteron gives S z = 0, what is the conditional probability for the next 
deuteron to pass through the Stern-Gerlach spectrometer also to be measured to have S z = 0? The answer this time 
is no longer 1/3. The spins of all the deuterons are aligned along the same direction n, and that the first deuteron is 
measured to have S z = suggests n is more probable to be more or less aligned along z than otherwise. As a result, 
the conditional probability is no longer 1/3, but can be easily shown to be 3/5, which is exactly the predicted value 
for (R)i in Eq. (^p|). The situation for a single domain of D%C is analogous, with isospin aligned pions instead of spin 
aligned deuterons. By construction, the pions in a D%C domain are isospin aligned, and by the same analysis, we 
have shown that the knowledge of the "first pion" being neutral can dramatically modify the conditional probability 
distribution of R. 

One can also consider the conditional probability distribution of R in the case that the "first pion" is charged. Note 
that 

fo(R) — |(/(-R|lst pion is a ir + ) + /(i?|lst pion is a 7r_) + /(i?|lst pion is a 7To)), (11) 

and hence, since fi{R) = /(i?|lst pion is a n ), 

f(R) = /(i?|lst pion is charged) = §/ (#) - \h{R) = |(1 - R)R- 1 / 2 . (12) 

The expectation value of R, given that the "first pion" is charged, can be easily shown to be 1/5. As a consistency 
check, one can calculate (R) , the expectation value of R regardless of the species of the "first pion" . Since the "first 
pion" is twice as likely to be charged as to be neutral, 

(#>o = i(f + 2xi) = i, (13) 

agreeing with Eq. (Q). 

Lastly, we will study the conditional probability distribution of R given that the k pions with the lowest pr, which 
will be hereafter referred to as the "first k pions" , are all neutral. It is straightforward to show that in this case 



and 



dp = 2— — = l Rk - R 1/2dR (u) 

J±sm8cos 2k 9de Lr* . R-^dR 



f k {R) = f(R\lst k pions are all neutral) = R k f {R) j R k fo{R)dR = (k + \)R k ~ x / 2 . 



(15) 



The distributions f2(R) and /a(i?) are plotted in Fig. lc and d, respectively. One can see that as k increases, the 
distribution is more and more skewed towards the high end. As a result, the expectation value of R increases with k. 
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(R) k = [ Rf k (R)dR=(2k+l)/{2k + 3). 
Jo 



(16) 



It is useful to define Q as the ratio of the number of 7r+ to the number of total pions emitted. By symmetry it is also 
the ratio of the number of 7r_ to the number of total pions emitted, and since R + 2Q = 1, 

(Q) fc = l/(2fc + 3). (17) 

From the above analysis, the prescription to enhance the collective signal is quite clear. One should make successive 
cuts on the data sample on the condition that the k pions with the lowest pt are all neutral, and measure (R)k after 
each cut to see if it increases as predicted in Eq. (|l6|). This result, however, depends on the assumption that we have 
only a single domain of D^C without any contamination due to incoherent pion emissions. Since this assumption is 
unrealistic for heavy ion collision experiments, the scenario we studied in this section is only an idealized situation. 
In the next section, we will discuss more realistic scenarios. 

III. THE EFFECTS OF MULTI-DOMAIN FORMATION AND INCOHERENT EMISSIONS 

The scenario considered in the last section is highly unrealistic in at least two ways. First, as discussed in the 
introduction, single domain D%C formation is highly unlikely. For a realistic treatment one must study D^C formation 
with more than one domain, each pointing in a different direction in the isospace. Moreover, we have neglected the 
effect of incoherently emitted pions, which have very important effects. If the neutral "first pion" is incoherently 
emitted, the R distribution of the remaining pions is described by fo(R), instead of fi{R) when the "first pion" comes 
from the D^C. In this section, we will incorporate these two effects and see how the predictions above are modified. 

We will study the expectation value of R, or equivalently the expectation of Q, for a situation described by the 
following parameters. The coherent fraction \ is the fraction of pions which originate from D^C domains, so that 
when x = 1, all pions are coherently emitted, and when \ — 0, all pions are incoherently emitted. We will consider the 
case where there are N domains, all containing an equal number of pions^J, which will be assumed to be large. Each 
domain is described by an isosinglet density matrix, but the isospins of pions in different domains are uncorrelated. 
Now the question is: if the "first k pions" in this channel are all neutral, what are the expectation values of R and Q 
among the rest of the pions? 

The answer turns out to be the following expression: 

§ - A. (18a) 



j=0 J 



V = X/N. (18c) 

Each term in this formula has a simple interpretation: 

• The expectation value (R) is always 1/3 for the incoherently emitted pions. Only the pions coming from the 
domains are affected by isospin alignment; hence the outstanding factor of \- 

• Each coherently emitted pion comes from one of the domains, which will be called domain X. How many of the 
"first k pions" also come from domain X? The probability for each pion coming from domain X is p = x/iV, and the 
probability that j of the "first k pions" coming from domain X is Pj = ( .) p> (1 — p) k ~i . 

• Given that j of the "first k pions" is coming from domain X, the conditional expectation value of Q decreases 
from 1/3 to l/(2j + 3), while the conditional expectation value of R increases by twice the above quantity. 



3 This assumption of all domains having the same number of pions is unrealistic but is made for illustrative purposes. The 
effects of unequal domain sizes will be briefly discussed below. 



(R) = i + 2A, (Q) = 
The shift A is given by 



where 



p j = ( k )p i 0-p) k - i , 
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In passing, we note that A can also be expressed as an integral or the hypergeometric function 2F1' 



A = x(~-^=J o dzz 2 (z 2 + l-p) k ) 



f l cos 2fc + 3 Q f e \ a sin 2 ?? x 2 x/N 

XU rinr- / d # I~5feZJTo . tan 6 = 



^3 sin 3 9 7 cos 2fc + 4 ??y' ~ 1-X/Af 

= 5x(i-(i-P) fc a F l[ |,-*,|,^]). (19) 

In Fig. 2, we have made contour plots of A = 1/30 (such that (R) = 0.4 and (Q) — 0.3), for k — 1, ... ,5 in the 
(x, parameter space. The horizontal axis is the coherent fraction \ while the vertical axis is 1/N where N is 

the number of domains. Both \ and 1/7V range from to 1. Thus, for example, with 3 domains and x = 0-6, in order 
to have A > 1/30 we must have k > 3. 

Equations (Eq) illustrate the main results of this paper. One can see that, without any D%C formation, x = 
(corresponding to the left edge of Fig. 2), A vanishes, and (R) = (Q) — 1/3 as expected. The bottom edge of the plot 
corresponds to N — > 00 and also gives A = for any finite value of k. The shift A is largest for a single domain of 
D^C without any noise due to incoherently emitted pions, i.e., when y = N = 1 (the top right corner of the contour 
plots), giving A = 1/3 — l/(2fc + 3) and reproducing Eqs. ( |l6| ) and (|l7|). For fixed values of (x, N), A increases with 
k, accounting for the spreading of the parameter space with A > 0.4 as k increases from 1 to 5 in Fig. 2. 




FIG. 2. Contour plots of A = 1/30 (such that (R) = 0.4 and (Q) = 0.3) for different values of k. The horizontal axis is the 
coherent fraction x> while the vertical axis is 1/N where iV is the number of domains. Both \ and 1/N range from to 1. The 
curves, from top right to bottom left are for k = 1, 2, 3 4 and 5, respectively. The shift A is larger than 0.4 above the curves 
and smaller than 0.4 below the curves. 



One expects that when the number of D^C domains is large (N 3> 1) or when most of the pions are incoherently 
emitted (x <C 1), it will be difficult to observe clear signals of D%C formation. However, in such situations x/iV is 
small and A is dominated by the j ' = 1 term (the j ' = term always identically vanishes) and 
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Thus a large k may make up for a small coherent fraction x, or a large number of domains N, and enhance A, which 
describes the shift of (R) and (Q) from 1/3, to an experimentally measurable magnitude. From the form of Eq. (|2(i|), 
one expects this shift to be substantial whenever k ~ N/x 2 - However, even for a value of k as small as N/Ax 2 , 
A = 1/30 + 0(x 3 /N 2 ), which translates to (R) = 0.4 and (Q) = 0.3 — a substantial deviation from the incoherent 
case. This suggests one should make successive cuts for events where the k pions with lowest pr are all neutral, and 
study (R) after each cut. An increase of (R) with k would suggest that D^C domains are formed. 

Equation ( p0| ) appears to suggest that one can increase A to an arbitrarily large magnitude by choosing a sufficiently 
large value of k. Of course this is not true. Equation ( |2(i| ) is obtained as the leading term in a x/N expansion, but 
when k — ► oo, this expansion breaks down as terms of higher order in x/AT are enhanced by factors of Qj. We can 
easily see that 

A^ix, (#}-»! + §X. (Q)^i(l-x), fc -» oo with x and A fixed. (21) 

In other words, our signal enhancement scheme is fundamentally limited by the amount of noise due to incoherently 
emitted pions. When x is small, most of the pions are incoherently emitted, and for them, (R) is always around 1/3 
regardless of what cuts one makes. On the other hand, the large k limit of (R) does not depend on N, the number of 
D%C domains. Recall that we have several distinctive signatures, like the R distribution in Eq. ([!]) and the conditional 
expectation values for R described in the previous section, for a single domain of DxC, where all the pions in the 
D^C are isospin aligned. With multi-domain formation, where the pions in different domains may point to different 
directions in isospace, the effect of isospin alignment is greatly washed out. However, given that the "first k pions" 
are all neutral with fc 3> N, it is probabilistically extremely likely that each of the N domains is the origin of some 
of these "first k pions" . As a result, each of these N domains are well-aligned along the ttq direction, and hence also 
well-aligned with each other. As k — > oo, the N domains look more and more like a single big domain in the ttq 
direction, which is the case where the signature is the most dramatic. In short, with large k, our cuts are picking out 
the events where the signals are the strongest, and hence resulting in a large signal-to-noise ratio. 

In above we have assumed the sizes of all N domains are identical for illustrative purposes. A more realistic 
treatment would have N domains, each with different sizes pi, 1 < i < N, such that ^2 { Pi = X- (The size of a domain 
is defined to be the fraction of pions which originate from this particular domain.) Then Eqs. ( jig ) are generalized to 

^h-±*±0^-^^- 3 . (22) 

i=l j=0 VJ/ J 

In the weak signal limit, i.e., when all the p^'s are small, Eq. ( |20| ) gets modified to 

i—l 

where p — J^iPf nas the following nice interpretation: p is the average over all pions (both coherently and incoherently 
emitted) of the sizes of the originating domains, which is pi for a pion from domain i and zero for an incoherently 
emitted pion. For N domains of equal sizes, p = x 2 /A and Eq. ( pfj| ) is recovered. Again we see that A grows linearly 
with k in the weak signal limit. As k — > oo the shift A is again limited by the bound (pi]), which applies also for the 
cases of unequal domain sizes. 

IV. SUMMARY 

To recapitulate, we suggest the following experimental procedures: 

• Count the number of neutral and charged pions event by event from heavy ion collision experiments and measure 
their individual transverse momenta and rapidities. 

• Apply a low pt cut to suppress the noise due to uncorrelated pion emission. 

• Bin the events in different rapidity windows. 

• In each rapidity window, calculate the expectation value (R). 

• Make a cut to retain only events where the pion with the lowest px is neutral. 

• Calculate, in each rapidity window, the expectation value (R) for all remaining pions in all events which survive 
the cut. 

• Make another cut on the surviving events to retain only those where the pion with the second lowest pr is also 
neutral. 
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• Again, calculate in each rapidity window the expectation value (R) for all remaining pions in all events which 
survive the cuts. 

• Repeat the above prescription of making successive cuts to retain only events in which the pion with the next 
lowest pt is also neutral, and calculate (R) for each rapidity window after each cut. If we find (R) deviates from 1/3 
then we are seeing signatures from D^Cs. 

Note that this prescription requires reconstructions of pt's of individual pions, both charged and neutral. We have 
also presumed that the coherent fraction x an d the number of domains formed N are roughly the same for each event. 
(More specifically, the probability distributions for \ and N are narrow peaked.) 

By applying these successive cuts, we are retaining the events with D^C formation and most of the pions are 
well-aligned along the ttq direction. What is being cut are the events with D%C formation but most of the pions are 
well-aligned along the n x or n y directions, and the events where there are incoherent pions with very low pj>, which is 
the main source of noise to our signal. As a result, these successive cuts are substantially improving the signal-to-noise 
ratio, making it easier to observe D%C formation. On the other hand, just like any other cuts on data to suppress the 
noises, we are giving up on statistics. Moreover, for large k we are cutting on rare events so the loss in statistics can 
be severe. For the cases where the signal is weak (small coherent fraction x <C 1 or large number of domains N ^> 1) 
on each cut we are losing about two-thirds of the events. 

In conclusion, we have devised new cuts to enhance the signal in searches for D^C. These cuts retain only events 
where the k pions with lowest pr are all neutral. We have shown that, after these cuts, the fraction of neutral pions 
within the remaining sample is substantially larger if DxCs are formed in the heavy ion collision. 

Support of this research by the U.S. Department of Energy under grant DE-FG02-93ER-40762 is gratefully ac- 
knowledged. 
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